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Abstract. We present a many-body approach to the electronic and magnetic 
properties of the (multiband) Kondo-lattice model with ferromagnetic interband 
exchange. The coupling between itinerant conduction electrons and localized 
magnetic moments leads, on the one hand, to a distinct temperature-dependence 
of the electronic quasiparticle spectrum and, on the other hand, to magnetic 
properties, as e. g. the Curie temperature Tq or the magnon dispersion, which 
are strongly influenced by the band electron selfenergy and therewith in particular 
by the carrier density. We present results for the single-band Kondo-lattice model 
in terms of quasiparticle densities of states and quasiparticle band structures 
and demonstrate the density-dependence of the self-consistently derived Curie 
temperature. The transition from weak-coupling (RKKY) to strong-coupling 
(double exchange) behaviour is worked out. 

The multiband model is combined with a tight-binding-LMTO bandstructure 
calculation to describe real magnetic materials. As an example we present 
results for the archetypal ferromagnetic local-moment systems EuO and EuS. The 
proposed method avoids the double counting of relevant interactions and takes 
into account the correct symmetry of atomic orbitals. 



1. Introduction 

The Kondo-lattice model (KLM), being also denoted as s-f or s-d model or, in its 
strong-coupling regime, as double exchange model, is today surely one of the most 
prominent models in solid state theory, and that mainly because of its great variety of 
important applications to rather hot topics in the wide field of collective magnetism. 
It refers to magnetic materials which get their magnetic properties from a system of 
localized magnetic moments being indirectly coupled via an interband exchange J to 
itinerant conduction electrons. Many characteristic properties of such materials can 
be traced back to this interband exchange. 

There are a lot of important applications of the KLM, first of all to the classical 
systems of magnetic semiconductors such as EuO and EuS and magnetic metals 
such as Gd |2j (J > 0, weak to intermediate coupling). All these materials have strictly 
localized moments because of the half- filled 4/ shell of the rare earth ion (Eu 2+ , Gd 3+ ). 
The striking temperature-dependence of the unoccupied (!) conduction bandstates in 
the EuX Q] manifesting itself, e. g., in the well-known red shift effect PP, and the 
effective moment coupling in Gd via a conduction electron spin polarization (RKKY) 
can be understood only by the existence of the interband exchange J. 

Other candidates are the intensively investigated diluted magnetic semiconductors 
such as Gai-^Mn^As with probably negative J (strong coupling). The Mn 2+ ion 
creates simultaneously an S = | spin and an itinerant hole in the GaAs valence band 
H3- The latter takes care for an indirect coupling between the randomly distributed 
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magnetic moments that leads already for very small i to a ferromagnetic order with 
a maximum T c of 110 K for x = 0.053 0]. 

Another hot topic because of great potential for engineering applications are 
the colossal magnetoresistance (CMR) materials [5] such as Lai_ x (Ca,Sr) :c Mii03 
which show up a very rich magnetic phase diagram. For 0.2 < x < 0.4 
the originally antiferromagnetic and insulating parent compound LaMnOs changes 
into a ferromagnetic metal. This is ascribed to a homogeneous valence mixing 
(Mnfj^ Mn^ + ). The three 3dt 2g electrons in Mn 3+ are considered more or less 
localized forming an S — 1 moment while the additional 3d es electron is itinerant. 
Again, the local moment-itinerant electron exchange is made responsible for many 
typical features of the CMR materials (J > 0, strong coupling). 

We derive in the next section the multiband-Hamiltonian of the Kondo-lattice 
model which is combined in Section 4 with an "ab initio" bandstructure calculation 
to get the temperature-dependent electronic structure of the prototype ferromagnets 
EuO and EuS. To understand the physics of the KLM we present in Section 3 for the 
nondegenerate case an analysis of the magnetic and electronic model properties. 

2. Multiband Kondo-Lattice Model 

The system under consideration consists of quasi-free electrons in rather broad 
conduction bands (d bands) and localized electrons with extremely flat dispersions 
(f levels). While there is no contribution of the f electrons to the kinetic energy, the 
part of the d electrons reads 

Hd = Tij c tma C yma (1) 

ijmrha 

The hopping process from site Ri to site Rj may be accompanied by an orbital change 
(m — > m). T™ m are the respective hopping integrals. c- ma (cj ma ) is the creation 
(annihilation) operator for a Wannier electron at site Rj in the orbital m with spin a 
(a- =?,!)• 

When applying the model study to real local-moment materials, intended as our 
final goal, we have to require that the single-electron energies do not only account 
for the kinetic energy and the influence of the lattice potential, but also for all those 
interactions which are not explicitly covered by the model Hamiltonian. That means 
that the hopping integrals are to be taken from a proper "ab initio" band calculation. 
We exemplify the procedure in Sect. 4 for the special case of the ferromagnetic 
europium chalcogenides. In particular we show how to avoid the well-known double- 
counting problem of important interactions. 

The Coulomb interaction is restricted to intraatomic terms. Furthermore, we 
assume that only two subbands are involved in the scattering process. Then the 
interaction part can be written as |S] 

H c = H dd + H ff + H df (2) 

In an obvious manner the Coulomb interaction may be split into three different parts 
depending on whether both interacting particles stem from a conduction band H dd , or 
both from a flat band Hff, or one from a flat band and the other from a conduction 
band Hdf- The first term refers to electron correlations in the conduction bands, not 
explicitly taken into account by the KLM. As mentioned, our procedure for describing 
real materials accounts for such correlations by a proper renormalization of the single 
particle energies. 
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The term Hff describes interactions between electrons from the flat bands which 
interest us only with respect to the fact that they form permanent magnetic moments 
(spins): 

S i = J2^f (3) 

/ 

<r f is the spin operator of an electron in subband f. If it is necessary to consider 
a (super)exchange interaction between the local spins then Hff is chosen to be a 
Heisenberg Hamiltonian, possibly with a symmetry- breaking single-ion anisotropy D: 

H ff ^\/ (/ S,.S, ny^lS;r (4) 

ij i 

The third term in Hdf refers to the interaction between a localized and an 
itinerant electron. Neglecting unimportant spin-independent contributions it can be 
written as an intraatomic exchange, i. e. a local interaction between the conduction 
electron spin tTj m and the local moment spin Sj |5] : 

H d f = -J y^o-jm • Sj (5) 

jm 

Here m denotes the conduction electron orbital. J is the exchange coupling constant 
being assumed to be identical for all df pairs. Using second quantization for the 
itinerant electron spin (tij mcr = c~j mcr Cj mcr ), this interaction reads: 

H df = -\j Yj ( S j ( n J' m T " n omi) + S^cf ml c jml + Sjc+ m] c 3m ^j (6) 

jm 

We see that the first term describes an Ising-like interaction of the two spin operators 
while the two others provide spin exchange processes between localized moment and 
itinerant electron. In particular the two latter terms are responsible for many typical 
properties of the Kondo-lattice materials. Spin exchange may happen in three different 
elementary processes: magnon emission by an itinerant j electron, magnon absorption 
by an f electron and formation of a quasiparticle (magnetic polaron). The latter can be 
understood as a propagating electron dressed by a virtual cloud of repeatedly emitted 
and reabsorbed magnons corresponding to a depolarization of the immediate localized 
spin neighbourhood. 

The total Hamiltonian of the multiband-Kondo lattice is composed of , and 

m 

H = H d + Hff + H df (7) 

An important model parameter is of course the effective coupling constant where W 
is the Bloch-bandwidth. Especially the sign of J is decisive. Other parameters are the 
local spin S, the lattice structure, and above all the band occupation n = X)mo-( n mo-)- 
In case of a nondegenerate band n is a number in between and 2. 

3. Magnetic and Electronic Model Properties 

The many-body problem of the KLM is not exactly solvable for the general case even 
for the nondegenerate case, for which we can remove the orbital dependence and which 
shall exclusively be tackled in this section: 

H = J2 T v c tcj« ~jJ2 S r °i = H o + H! (8) 

ijo 3 
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However, there exists a non-trivial, very illustrative limiting case which is rigorously 
tractable and nevertheless exhibits all the above-mentioned elementary excitation 
processes [71IS1IH1E1- It refers to a single electron in an otherwise empty conduction 
band coupled to a ferromagnetically saturated moment system, e. g. EuO at T — 0. 
In this case the f spectrum is simple, because the f electron cannot exchange its spin 
with the parallel aligned spin system. Only the Ising-type interaction in J^J and © 
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Figure 1. l-spectral density as function of energy for several symmetry points 
in the first Brillouin zone and for different exchange couplings J. Parameters: 
S = §, W = leV, sc lattice. 



takes care for a rigid shift of the total spectrum by —^JS. The spectral density is a 
((-function at e(k) — \ JS where e is the free Bloch energy, the Fourier transform of 
the hopping integral IJy. Real correlation effects appear, however, in the J, spectrum. 
Fig. n shows the energy dependence of the |-spectral density for some symmetry 
points. For weak couplings the spectral density consists of a single pronounced peak. 
The finite width points to a finite quasiparticle lifetime due to first spinflip processes, 
but the quasiparticle bandstructure, being read off from the peak positions, is still 
very similar in shape to the Bloch dispersion. 

This changes drastically already for rather moderate effective exchange couplings 
JS/W . One observes in certain parts of the Brillouin zone, for a strongly coupled 
system even in the whole Brillouin zone, that the excitation energy splits into two 
parts. The sharp high-energy peak belongs to the magnetic polaron while the broader 
low-energy part consists of scattering states due to magnon emission. As long as the 
polaron peak is above the scattering spectrum the quasiparticle has even an infinite 
lifetime. The scattering spectrum is in general rather broad because the emitted 
magnon can carry away any wave-vector from the first Brillouin zone. Because of the 
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concomitant spinflip, magnon emission can happen only if there are f states within 
reach. Therefore, the scattering part extends just over that energy region where p-t ^ 
(p a : quasiparticle density of states (Q-DOS)). 

For the general case (finite temperature, finite band occupation) the many-body 
problem of the KLM cannot be solved exactly. Approximations must be tolerated. 
The central quantity is the selfenergy, for which we use the Green function procedure 
developed in 11 . This non-perturbational theory can be considered a "moment 
conserving decoupling approximation" (MCDA) which interpolates between exact 
limiting cases as, e. g., the above discussed example. In more detail, the method 
decouples the hierarchy of equations of motion by expressing certain "higher" Green 
functions as linear combination of "lower" functions, already involved in the procedure. 
The ansatz is motivated by exact special cases as the atomic limit, ferromagnetic 
saturation, local spin S = 1/2, full or empty band, and so on. Free parameters in 
the ansatz are eventually fitted to rigorous spectral moments. As other methods, too, 
e. g. the interpolating selfenergy approaches in refs.^J^li it leads to the following 
structure of the selfenergy: 

£k CT (£) = -\jZo (S z ) + J 2 D ka (E) (9) 

Restriction to the first term, only, yields the mean-field approach to the KLM, which 
is correct for sufficiently weak couplings J. It is due to the Ising-part in Eq. (J5J 
{z a = Scri — S a i). Without the second part of the selfenergy it would give rise 
to a spin-polarized splitting of the conduction band. The term D^ a {E) is more 
complicated being predominantly determined by the spin exchange processes. It is 
a complicated functional of the selfenergy itself, i. e. © is an implicit equation for 
Sker(-E') and not at all an analytical solution. D^^E) depends, further on, on mixed 

spin correlations such as {S* n ia ) Js^ cf^c,^ , built up by combinations of localized- 

spin and itinerant-electron operators. Fortunately, all these mixed correlations can 
rigorously be expressed via the spectral theorem by any of the Green functions involved 
in the hierarchy of the MCDA. 

However, there are also pure local-moment correlations of the form (S 1 ?), (5' I ± S'? = ), 

(Sf) 3 ^}, which also have to be expressed by the selfenergy 0. For this purpose 

we use the modified RKKY theory of |14| which exploits a mapping of the interband 
exchange (JSJ) to an effective Heisenberg model, 

H f = -J^J ij S i -S j (10) 

ij 

by averaging out the conduction electron degrees of freedom: 

- JJ2 Si • <tj — ■ W (C) — > H f ( n ) 

j j 

I + \ (c) 

In the last analysis this means to determine the expectation value (c£ + c kcr > ) . We 
use again a Green-function procedure Jl] which eventually leads to effective exchange 
integrals in l|lUfl : 

t2 r + oo 

4 = T^72E eiq ' (R!_Rj) / dEUE)lm[(E~e(k) +l 0+) 

kqa J -°° 

x (E-eik + qj-^+^E))}- 1 (12) 
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f-(E) denotes the Fermi function. The effective exchange integrals are decisively 
influenced by the conduction electron selfenergy E CT which brings a distinct band 
occupation and temperature dependence to the J^-. Neglecting H a leads to the 
"conventional" RKKY formula with J.^ ex; J 2 as a result of second order perturbation 
theory. Via E CT higher order terms of the electron spin polarization enter the modified 
RKKY being therefore not restricted to weak couplings, only. Note that at the right- 
hand side of Eq. H12(l the product of two propagators appear. From the equation of 
motion method in it follows directly that one of them is dressed by a selfenergy, 
the other not. 

To get from the effective operator <|10|) the magnetic properties of the KLM 
we apply the standard Tyablikov-approximation which is known to yield convincing 
results in the low as well as in the high temperature region All the above- 

mentioned local-moment correlations are then expressed by the electronic selfenergy. 
We therefore end up with a closed system of equations that can be solved self- 
consistently for all entities of interest. 

Fig. [21 shows the temperature dependence of the quasiparticle bandstructure (Q- 




Figure 2. Quasiparticle bandstructure as a function of wave- vector (left column) 
and quasiparticle density of states as a function of energy (right column) for four 
different temperatures calculated within the MCDA lTTl . Parameters: J = 0.2 eV, 
W = 1 eV, n = 0.2, S = |, sc lattice. 

BS), derived as density plot from the spectral density, and the quasiparticle density 
of states (Q-DOS) p a {E) for a typical parameter set from the moderate coupling 
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regime (J — 0.2 eV, W — 1 eV, n = 0.2, S = |, sc lattice). In this case the self- 
consistently calculated Curie temperature Tq amounts to 238 K. At T — 37 K the 
local-moment magnetization is 3.4 and therewith very close to saturation. The theory 
reproduces, in spite of the finite carrier concentration, the same features as for the 
exact (n = 0, T = 0)-case exhibited in Fig. ^ The f dispersion as well as the QDOS 
pij are practically identical to the respective "free" Bloch functions. The J, spectrum, 
however, is more complicated. The scattering states take away a substantial part of 
the spectral weight near the T point while near the R point the polaron state clearly 
dominates. With increasing temperature (decreasing (S z )) a finite magnon density 
appears allowing for scattering states in the f spectrum, too, because of magnon 
absorption and simultaneous spinflip by the | electron. Furthermore, polaron and 
scattering states separate, where, surprisingly, the rather broad scattering spectrum 
is in wide wave-vector regions bunched to a prominent peak. At T = Tq the spin 
asymmetry is removed but there remains a correlation-caused splitting of the spectrum 
into two branches and two subbands, respectively, due to the interband exchange 
coupling J. Because of the possibility of mutual spin exchange and pi are occupying 
for finite temperatures always the same energy regions. 

The key-quantity of ferromagnetism is the Curie temperature Tq. A finite Tq 
comes out as a consequence of an indirect coupling between the local moments, 
mediated by a polarization of the conduction electron spins. Therefore, Tq 
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Figure 3. Curie temperature as a function of (a) the band occupation n for 
various J, (b) the interband exchange coupling J for various n, calculated by use 
of the modified RKKY. 



exhibits a strong carrier-concentration dependence which is plotted in Fig. [3J The 
effective exchange integrals have been taken into account up to 25th neighbours. 
Ferromagnetism appears for low electron (hole) concentrations, while being excluded 
in the region around half- filling (n = 1). Around n = 1 we expect an antiferromagnetic 
phase, however, not explicitly calculated here. A similar n-dependence of Tq has been 
found in Ref. |17| . It is interesting to compare the results with those of the conventional 
RKKY, given by = in Eq. Ijl2|l . A corresponding curve is inserted in Fig. The 
maximal Tq values are higher, but ferromagnetism exists only in a very narrow region 
of low electron (hole) concentrations where this region turns out to be independent of 
J. 

Two features dominate the J-dependence of the Curie temperature plotted in 
part (b) of Fig. The first is the appearance of a critical J for band occupations 
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for which the conventional RKKY does not allow ferromagnetism (n > 0.13). This 
is in accordance with the fact that for J — > the modified RKKY reproduces the 
conventional RKKY. The second fact is that Tq runs into a saturation for strong 
couplings J, where the saturation value depends on the band occupation n. These 
results are far beyond conventional RKKY which can work of course only in the weak 
coupling regime. 

It is interesting to observe (see Fig. that in the weak coupling RKKY region 
Tc scales with the effective coupling constant and in the strong coupling double 




Figure 4. Curie temperature as a function of (a) the interband exchange J 
for various W and n = 0.1, (b) the Bloch bandwidth W for various J and n, 
calculated by use of the modified RKKY. 



exchange region with the kinetic energy, i. e. oc W. In part (b) of Fig. 0] the change 
from RKKY behaviour (small ^) to double exchange behaviour (large -d?) is clearly 
to be seen. 

Fig. shows a typical example of a spin-wave dispersion, which obtains via the 
effective exchange integrals lfT^|) a distinct temperature dependence. Upon heating the 
dipersion relation uniformly softens, disappearing above Tc. This agrees qualitatively 
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Figure 5. Spin-wave dispersion of the ferromagnetic Kondo-lattice model as a 
function of the wave-vector for different temperatures T. Parameters: J = 0.2 eV, 
n = 0.2, 5 = 7/2, sc lattice, and W = 1 eV. The sclf-consistently calculated Curie 
temperature is Tq = 232 K. 
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with neutron scattering data on manganites |18| . 

4. Application to ferromagnetic EuO and EuS 

In the last step we want to demonstrate how the model analysis can be used to describe 
real magnetic materials. We exemplify the method with respect to the classical local- 
moment ferromagnets EuO and EuS. These are ferromagnetic semiconductors which 
crystallize in the fcc-rocksalt structure with lattice constants a = 5.142 A(EuO) and 
a = 5.95 A(EuS), respectively. The physics of these materials is mainly determined 
by localized and half-filled 4f levels, on the one hand, and, on the other hand, by 
extended, at low temperature empty, 5d band states. In accordance to Hund's rule, 
the Eu 2+ -4f shell produces a maximum spin moment S = |. 

To investigate real materials such as EuO and EuS one has to exploit the 
multiband version of the KLM developed in Section 1. As already explained, the 
hopping integrals T^ m in have not only to incorporate the kinetic energy of 
the band electrons plus the influence of the periodic lattice potential but also as 
realistically as possible all those interactions which are not directly covered by the 
model Hamiltonian. For this purpose we performed a bandstructure calculation 
based on density functional theory, using the Andersen scheme EH1 of a tight- 
binding linear muffin-tin orbital (TB-LMTO) ansatz. LDA-typical difficulties arise 
with the strongly localized character of the 4f levels. A normal LDA calculation for 
EuO produces a metal with the 4f levels lying well within the conduction band. To 
circumvent the problem we considered the 4f electrons as core electrons, since our main 
interest is focused on the temperature-reaction of the (empty) conduction bands. In 
our study the 4f levels appear only as localized spins (moments) in the sense of H/f 
in Eq. QJ. Fig. shows the so-derived spin-dependent bandstructure of EuS without 
the 4f levels. The low-energy part belongs to S-3p states, we are not interested in 
here, because only fd excitations shall be in the focus. The conduction band region 
is dominated by Eu-5d states, to which we have restricted the investigation by our 
combined many-particle, first principles procedure. For comparison we have also 
performed an LDA+U calculation which indeed is able to put the respective bands 
into better positions. On the other hand, the results, which we are interested in, do 
not differ too strongly from those of LDA with the 4f electrons in the core. So we 
have chosen the much simpler LDA calculation, since our study mainly aims at overall 
correlation and temperature effects. The extreme details of the band structure are 
not so important. 

The many-body treatment of the multiband KLM follows exactly the same line 
as for the single-band case, no need to create new approaches. However, we have to 
fix the only parameter of our procedure, the interband exchange J. It is generally 
assumed that an LDA treatment of ferromagnetism is quite compatible with a mean- 
field picture, so that the T = 0-shift of the j and J, dispersions in Fig. should amount 
to A = JS. However, it can be seen that the assumption of a rigid shift is too simple. 
A certain energy dependence of the exchange splitting is found by LDA, too. We 
decided to use the splitting of the 5d-band centre of gravity to fix the J for EuS at 
0.23 eV PS]- Averaging over prominent features in the Q-DOS of EuO yields a similar 
value of J = 0.25 eV [E]. 

Another important issue concerns the renormalized single-particle input, because 
a double counting of the decisive interband exchange, once by the LDA input and 
once more explicitly in the model-Hamiltonian might produce rather misleading 
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Figure 6. Spin-dependent (solid lines down-spin, broken lines up-spin) band 
structure of bulk EuS calculated within a TB-LMTO scheme with the 4f levels 
treated as core states. The energy zero coincides with the Fermi energy. 



results. The most direct solution of this problem would be to switch off the interband 
exchange Hjf in the LDA code, which turns out to be impossible. However we can 
elegantly exploit the exact limiting case, discussed at the beginning of Section 3. For 
zero temperature and zero band occupation the | spectrum of the single- as well as 
multiband-KLM is only rigidly shifted towards lower energies by the amount of \ JS. 
Thus we can use without any manipulation the f dispersions of the LDA calculation 
(T = 0!) for the empty 5d bands as input for the single-electron Hamiltonian (JIJ. There 
is no need to switch off Hjf because in this special case it leads only to an unimportant 
rigid shift. Note that this holds only for the f part, the [ spectrum is even at T = 
strongly influenced by correlation effects (see Fig. QJ. Therewith it is guaranteed, on 
the one side, that all those interactions, which are not directly covered by the KLM, 
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are implicitly taken into account by the LDA-renormalized single-particle energies. On 
the other side, a double counting of any decisive interaction, which usually occurs when 
combining first-principles and many-body model calculations, is definitely avoided. 

Because of the empty conduction bands a self-consistent justification of the 
EuO (EuS) ferromagnetism via the modified RKKY-treatment of Section 3 is not 
possible. The magnetic properties of the insulators EuO and EuS are exclusively 
due to the partial operator Hff (@J. Its exchange integrals between nearest (Ji) and 
next-nearest {J 2) neighbours have been derived by a neutron scattering spin wave 
analysis |H]: Ji/fc B = 0.625(0.221) K, J 2 = 0.125(-0.100) K for EuO (EuS). The 
single-ion anisotropy constant D does not play a decisive role as long as the bulk 
bandstructure and its temperature dependence is aimed at. However, when treating 
systems of lower dimensionality (films, surfaces) as done for EuO in Refs. [2U El, 
then a finite D is necessary to overcome the Mermin- Wagner theorem 22 for getting 
a collective magnetic order of the moment system. D has been chosen in such a 
way that the experimental Tq value is correctly reproduced (see Fig. [3] in Ref. |15|. 
The temperature-dependence of the 4f moment system leads to a remarkable induced 
temperature behaviour of the (unoccupied!) 5d spectrum. The evaluation of the 
effective Heisenberg model (0J is done by a Tyablikov decoupling of the equation of 
motion of a poperly chosen spin Green function. Details are given in |14l 111) . 

Let us now select some typical results found for the archetypal ferromagnetic 
semiconductors EuO and EuS. We have no intention to be complete in this overview, 
rather to demonstrate the possibilities of our combined many-body/afe initio method. 



<S Z >/S = 1, 0.75, 0.5, 0.25, 




Energy (eV) 

Figure 7. Temperature-dependent densities of states of the Eu-5d bands of bulk 
EuO (J = 0.25 eV). Solid lines: f-spectrum; broken lines: J-spectrum. Thick line: 
spectra for T = T c ({S z ) = 0). The outermost curves belong to T = ((S z > =5). 
With increasing temperature the | and J. curves approach each other. 



Fig. and |Hl show the quasiparticle density of states (Q-DOS) for EuO and EuS, 
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Figure 8. The same as in Fig. but for bulk EuS (J = 0.23 eV). 



respectively, for five different 4f magnetizations, i. e. five different temperatures. 
Not surprisingly, the Q-DOS of EuO (Fig. 0) and EuS (Fig. [BJ are qualitatively 
rather similar. Since we have taken into account the full bandstructure of the Eu- 
5d conduction bands, the symmetry of the different 5d orbitals is preserved. The 5d 
bands can therefore be decomposed into ti g and e g subbands, where the ti g bands 
are substantially broader (~ 7 eV (EuS), ~ 9 eV (EuO)) than the e g bands (~ 4 eV 
(EuS), ~ 6 eV (EuO)). A remarkable temperature-dependence shows up for both 
ferromagnets, manifesting itself, e. g., in a shift of the lower edge of the f Q-DOS to 
lower energies upon cooling from T = Tq to T = 0. This explains the famous red shift 
of the optical absorption edge for the electronic 4/ — 5d t2g transition, in the meantime 
observed for all ferromagnetic local-moment systems [23] • The temperature behaviour 
do not comply with the simple Stoner picture of an energy-independent induced 
exchange splitting of the spectra. The reason for the observed more complicated 
behaviour as a function of temperature is that in our theory correlation is treated in 
a way distinctly beyond mean-field. 

While the Q-DOS refers to the spin resolved but angle averaged, direct or inverse 
photoemission experiment, the k dependent spectral density is the angle resolved 
counterpart. Prominent peaks of the spectral density as a function of energy determine 
the E(k) quasiparticle bandstructure (Q-BS). As an example, Fig. |5] represents the 
EuO-spectral density as a density plot for some symmetry directions and three different 
temperatures. The degree of blackening is a measure of the magnitude of the spectral 
density. For T = ((S z ) = S) the t spectral density agrees, except for a constant 
shift, with the LDA result, being therefore a <5-function at the excitation energy, 
corresponding to a quasiparticle with infinite lifetime. The J, spectral density, however, 
exhibits already at T = a broadening of the dispersion curves, most notably around 
the T-point, which indicates finite quasiparticle lifetimes due to correlation effects. For 
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Figure 9. Spin-dependent quasiparticle band structure of the 5d-bands of bulk 
EuO (J = 0.25 eV) for three different 4f magnetizations {S z ). 



intermediate temperatures ((S z ) = 0.5 S) a broadening sets in also for the | dispersions 
because the f electron, too, can now exchange its spin with the no longer saturated 
local moment system. In the J, spectrum correlation effects increase leading to further 
quasiparticle damping. Simultaneously, both spin spectra are shifted towards each 
other reducing therewith the effective exchange splitting. At T = Tq {{S z ) = 0) f and 
J. dispersions conincide. Very similar behaviour is found for EuS |15| . 

As mentioned, the Q-BS has been derived from the spectral density. To give 
an insight into the temperature behaviour of the spectral density we have plotted in 
Fig. 1101 for two k-points from the first Brillouin zone (r, W) the energy dependence 
of the spectral density for the same three 4f magnetizations as in Fig. |5J but now 
for EuS. For the T-point (Fig. (a)) we have to expect two structures according to 
the twofold degenerate (e s ) and threefold degenerate fog) dispersions. Well defined 
quasiparticle peaks appear, spin split below Tq. The spin splitting, induced via the 
interband coupling to the magnetically active 4f system, collapses for T — > Tq. One 
speaks of a Stoner-type behaviour. Quasiparticle damping, that scales with the peak 
width, obviously increases with temperature. At the W-point (Fig. [23 (b)) four sharp 
peaks show up in the | spectrum, and, though already strongly damped, the same 
peak sequence appears in the J. part, too (compare with the Q-BS in Fig. |^. The 
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reduction of the spin splitting with increasing temperature leads to a strong overlap 
of the two upper peaks, which are then no longer distinguishable. Altogether, the 





Figure 10. (a) Spin-dcpcndcnt spectral density S^ct for k from the T-point (5d 
states) of bulk EuS (J = 0.23 eV) as a function of the energy for three different 4f 
magnetizations: (S z )/S = 1.0,0.5,0. Thick line for T = Tq. Solid lines: up-spin: 
broken lines: down-spin. Curves with maximum spin splitting of the peaks belong 
to T = 0. (b) The same as in (a) but for the W-point. 



5d spectral densities of the ferromagnetic semiconductors EuS and EuO (not shown 
here) exhibit drastic temperature-dependencies, what concerns the positions and the 
widths of the quasiparticle peaks. The driving force is the interband (4f-5d) exchange 
interaction, the investigation of which was the main aim of the presented study. 

5. Summary 

We have demonstrated how the (multiband) Kondo-lattice model can be used for 
a realistic description of the strikingly temperature-dependent electronic excitation 
spectrum of local-moment ferromagnets. For this purpose we have combined a many- 
body evaluation of the KLM with a first principles bandstructure calculation (TB- 
LMTO) for the ferromagnetic semiconductors EuO and EuS. 

In the first step we inspected an exactly solvable, non-trivial limiting case of the 
KLM to learn about the fundamental elementary excitations. In addition, this limiting 
case provides a weighty testing criterion for unavoidable approximations in the case 
of the not exactly solvable general case. 

By use of a self-consistent approach we found out the electronic as well as magnetic 
properties of the KLM as functions of decisive model parameters such as band 
occupation and interband exchange coupling. By inspecting the Curie temperature 
we could demonstrate the transition from a weak-coupling RKKY behaviour to a 
strong-coupling double exchange scenario. It appears remarkable that both features 
are simultaneously involved in our theory. 

Using a straightforward multiband version of the KLM with identical theoretical 
approaches we have worked out the influence of the localized magnetic 4f moments of 
the Eu 2+ ion on the 5d conduction bands in EuO and EuS, respectively. The normally 
empty bands exhibit a distinct and unconventional {"not Stoner-like") temperature- 
dependent exchange splitting induced by the interband exchange coupling to the 4f 
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levels. As a special detail we observe the famous red shift effect in EuO as well as 
EuS. 
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